Duality, Quantum Skyrmions and the Stability of an SO (3) 2D Quantum Spin-Glass 



o 
o 

(N 



(N 



CO 



i 

c 

o 
o 



> 
o 

o 

OS 

o 



C.M.S. da Conceigao 1 and E.C.Marino 2 

1 Departamento de Fisica Tedrica, Universidade do Estado do Rio de Janeiro, Rio de Janeiro, RJ 20550-013, Brazil and 

2 Instituto de Fisica, Universidade Federal do Rio de Janeiro, 
Cx.P. 68528, Rio de Janeiro, RJ 21941-972, Brazil 
(Dated: July 27, 2009) 

Quantum topological excitations (skyrmions) are analyzed from the point of view of their duality 
to spin excitations in the different phases of a disordered two-dimensional, short-range interacting, 
SO(3) quantum magnetic system of Heisenberg type. The phase diagram displays all the phases, 
which are allowed by the duality relation. We study the large distance behavior of the two-point 
correlation function of quantum skyrmions in each of these phases and, out of this, extract infor- 
mation about the energy spectrum and non-triviality of these excitations. The skyrmion correlators 
present a power-law decay in the spin-glass(SG)-phase, indicating that these quantum topological 
excitations are gapless but nontrivial in this phase. The SG phase is dual to the AF phase, in the 
sense that topological and spin excitations are respectively gapless in each of them. The Berezinskii- 
Kosterlitz-Thouless mechanism guarantees the survival of the SG phase at T / 0, whereas the AF 
phase is washed out to T = by the quantum fluctuations. Our results suggest a new, more sym- 
metric way of characterizing a SG-phase: one for which both the order and disorder parameters 
vanish, namely (a) = 0, (fi) — 0, where a is the spin and fi is the topological excitation operators. 

PACS numbers: 75.50.Lk 



I. INTRODUCTION 

Given a physical system, the spectrum of its possible 
excitations can, in principle, be roughly divided into two 
groups: hamiltonian and topological excitations. The 
first ones correspond to states, which are created out of 
the ground state by the action of operators that appear 
explicitly in the hamiltonian. These excitations carry 
in general quantum numbers corresponding to physical 
quantities, such as charge for instance, which are con- 
served as a consequence of some continuous symmetry of 
the hamiltonian. Topological excitations, conversely bear 
quantum numbers that correspond to quantities whose 
conservation derives from a nontrivial topology of the 
space of classical configurations and whose continuity 
equation is satisfied as an identity rather then being the 
result of some continuous symmetry in the system. 

In this work, we investigate the correlation functions 
of quantum skyrmions, which are the topological excita- 
tions that may occur in two-dimensional magnetic sys- 
tems described by a hamiltonian of the Heisenberg type, 
with SO (3) symmetry and nearest-neighbors- interactions 
on a square lattice. We consider the quenched disordered 
case, with a Gaussian random distribution of exchange 
couplings centered at an antiferromagnetic (AF) coupling 
J > 0, with variance A J. 

This disordered system has been studied by means 
of a mapping onto a generalized nonlinear sigma model 
(NLSM) , in which the original (staggered) spin is mapped 
onto the NLSM field n Q [2, H| and the index a = 1, ...,n 
corresponds to the different replicas which are required 
for a quenched average. Eventually we must take the 
limit n -> [11111. 

A remarkable feature of this system is that, whereas 
the quantum fluctuations completely "wash out" , the 



ordered Neel phase at T ■£ 0, in agreement with the 
Mermin- Wagner theorem [fa], the spin glass phase per- 
sists even at a finite temperature. Even though this does 
not violate the theorem, because there is no spatial long 
range order in the SG phase, it is somewhat intriguing 
to have a SG phase at a finite temperature in a quan- 
tum 2D system, specially if we consider a lot of evidence 
against the occurrence of a SG phase at a finite tempera- 
ture in 2D Ising systems [7[. We will see below, however, 
as a consequence of our study of quantum topological 
correlation functions, that this fact can be understood 
as a manifestation of the Berezinskii-Kosterlitz-Thouless 
(BKT) mechanism US]. Indeed, it is the BKT mecha- 
nism that allows the existence of a SG phase at a finite 
temperature in this 2D system. Moreover, the action 
of this mechanism is only possible in the case of SO (3) 
systems, where a vortex picture for the skyrmions does 
exist. This explains why the corresponding SG phase is 
not found in 2D Ising sy terns, where such a picture is 
absent. 



II. DUALITY AND QUANTUM TOPOLOGICAL 
EXCITATIONS 

The non-triviality of the topology of classical config- 
urations space, is ultimately responsible for the stabil- 
ity of classical topological excitations. At the quantum 
level, this nontrivial topology manifests as a degeneracy 
of the ground state. The existence of ground-state de- 
generacy, therefore, is the indication that, at a quantum- 
mechanical level, the system presents quantum topolog- 
ical excitations in its spectrum. These excitations, dif- 
ferently from the former, cannot be created by acting on 
the ground state with hamiltonian operators. A familiar 
example of topological excitations are magnetic vortices 
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in 2D, which occur in type II superconductors. The cor- 
responding topological "charge" would be the magnetic 
flux piercing the plane. 

Topological excitations, surprisingly, are related to the 
concept of order-disorder duality, which plays an impor- 
tant role in many areas of physics. This can be seen in 
the following way. Consider a system characterized by a 
dynamical variable a. It could be, for instance, the mag- 
netization on a ferromagnetic system or the staggered 
magnetization in an antiferromagnetic one. Suppose the 
hamiltonian H[a] possesses a symmetry, which at the 
quantum level is implemented by an unitary operator U 
such that, for g being the element of the symmetry group, 
UoW = ga, with [H, U] = 0. Since 



= 5 <o|t/W|o), 



(2.1) 



it follows that if (a) ^ 0, then necessarily U\0) = |0') ^ 
|0), that is, the ground-state will be degenerate. This, 
however, as remarked above, is a sign of non-triviality 
of the topological excitations. Now, assume the one- 
particle quantum topological excitation state is given by 
| Top) = /z|0), where /z is the operator that creates these 
excitations out of the ground state. Then, if these states 
are nontrivial, we must have (0|/i|0) = because this 
means |Top) is orthogonal to the ground state, as a gen- 
uine excited state must be. We conclude, therefore, that 
(a) ^ implies (/x) =0. 

Conversely if (0|/z|0) ^ 0, this would mean that |Top) 
is not orthogonal to the ground state and hence is ac- 
tually not a genuine excitation. This would imply the 
ground state should be unique, because otherwise topo- 
logical excitations should exist as nontrivial states. Now, 
if the ground state is unique, we have U\0) = |0) and 
(|2.ip would imply (a) =0. We therefore conclude that 
(/z) ^ implies (a) =0. 

The previous analysis shows that, if (a) measures the 
amount of order (an order parameter) then (/z) measures 
the amount of disorder, being naturally a disorder pa- 
rameter. We see that a duality relation exists between 
the topological excitations creation operator /z and the 
hamiltonian operator a. The physical and mathemati- 
cal properties of this duality relation are captured by the 
so-called dual algebra, satisfied by a and /x namely 



/"(x, *)cr(y, t) = g(y - x)cr(y, £)/j(x, t), 



(2.2) 



where, for each fixed x and y, g(x — y) is an element 
of the symmetry group. This relation is the basis for 
constructing the topological excitations creation operator 
zz. It implies, for instance, 



{*){») = 0, 



(2.3) 



hence we cannot have both (a) and (zz) nonvanishing. 
It also implies the spectrum is gapless whenever both 
(a) = and (/z) = [HI- According to (USD, for these 
kind of systems (with just one scalar order parameter) we 
can basically have only three phases. One with (a) ^ 



and (/z) = 0, another with (a) = and (zz) ^ and 
finally a phase with both (a) — and (zz) =0. 

From the large distance behavior of the quantum topo- 
logical correlators, we show that the three phases allowed 
by the dual algebra are realized in the disordered system 
considered below. 



III. THE CP 1 FORMULATION OF THE 
SHORT-RANGE AF HEISENBERG SPIN-GLASS 

A. The Nonlinear Sigma Model Formulation 

We review in this and in the next subsection the field- 
theoretical description of the disordered SO (3) quantum 
Heisenb erg- like system, which was studied with great de- 
tail in P, Q. It is described by the hamiltonian operator, 



(3.1) 



with nearest-neighbor interactions on a 2D square lattice 
of spacing a, having the random couplings associated 
with a Gaussian probability distribution with variance 
A J and centered in J > 0, namely, 



1 



V2t(AJ)' 



: exp 



(Jjj - Jf 

2(AJ) 2 



(3.2) 



We consider the case of quenched disorder, which is 
conveniently dealt with the replica method [H, |j| . In this 
case the average free-energy is given by 



F = -k B T lim - ([Z n ] 
n — >o n 



!)■ 



(3.3) 



where [Z n ] av is the disorder-averaged replicated partition 
function. 

With the help of spin coherent states \ flf (t)) such that 



(n?(T)|s?|n?(r)) = sn?(r) 



(3.4) 



where S is the spin quantum number, it was shown [l[ 
that in the continuum limit, the average replicated par- 
tition function corresponding to (|3.ip could be written 
as 



[Z n ]av 



J PnPQPAexp j- J drLj A [n a , Q a/3 , A Q ]| , 

(3.5) 

where the corresponding lagrangian density is a general- 
ized relativistic nonlinear sigma model (NLSM) 



i|Vn Q | 2 + ^|d T n«| 2 



iX a (\n° 



D 



dr' 



a/3, 



(3.6) 



-<(«f(T,T'K(T') 

Ps 
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where D = S 4 (AJ) 2 /a 2 (a: lattice parameter) and p s — 
S 2 J . In the above expression, summation on the replica 
indices a, ft = 1, n is understood. 

The field n Q = (er Q , 7? a ) is the continuum limit of the 
(staggered) spin f2 Q and satisfies the constraint n Q n" = 
p s , which is implemented by integration on A Q . 

Decomposing Q a @ into replica diagonal and off- 
diagonal parts, 

Q a J(r; r, r') = 6 ab [S^ X (r, r, r') + q^(r; r, r')] (3.7) 
where q a ^ = for a — (3, we get 



-J.A 



|Vn° 
3D 



2c 2 



\d T n a \ 2 +i\ a {\n° 



Ps) 



+ — I dr' [n X 2 {T,T')+ q ^{T,T') q ^{ry) 



-n a (r) X (r,r')n a (r') 

Ps 

— n Q (r)g Q/3 (r, r')n (r') 

Ps 



(3.8) 



This was the starting point for the CP 1 formulation, 
which was derived in [2| . In a previous work [l[ , we took 
a different path: from (13. 8j) . we integrated over the ff- 
field and thereby obtained an effective action for the re- 
maining fields. This allowed the determination of the 
average free-energy and, out of this, the phase diagram 
of the system. Identical results were found within the 
CP 1 framework [2|. 



B. The CP 1 Formulation 

Introducing the CP 1 field as usual, by the relation 
1 



n»(r) 



/p. 



[zr(r)a ijZ f(r)] 



where the zf field satisfies the constraint 



_a|2 i |„a|2 



1*1 I ' \ / -2 I 

and using the correspondence 



(3.9) 



(3.10) 



||Vn Q | 2 + ^\d T n*\ 2 ^2j2\D,zr\ 2 , (3.11) 



i=l 



(D^ — + iA^), which comprises a functional integra- 
tion over the auxiliary vector field A^, we obtain |2J 

[Z n ] av = J VzVz*VA^V X VqV\e- s , (3.12) 

where S [zf , zf *, A^, A, x(r, t'), q a0 (r, r')] is the action 
corresponding to the lagrangian density 

£j,A,c P1 = 2|iVf | 2 + ^ a (\zf\ 2 - p s ) 



+ ^ / dr' [n X 2 (r,r') + q^(r,r')q^(r,r')] 



y 'dr' {[ X (t,t>) 



l^r(T)| 2 |^(r')| 21 



- [zrzf(r)][ X (r,r')S^+ q ^(r 7 r')}[z^zf(r')}} , 

(3.13) 

where summation in i, j, a, [3 is understood. 

In we have determined the average free-energy, 
by expanding the fields around their stationary point in 
(|3.12[) and integrating the quadratic quantum fluctua- 
tions of the zf fields, namely 



S [ z f ' Z T* > A-n 7 X a , X ,q a ] — S [z" s , zf * , A s ^, m , Xs ,q s 



+ -J drdr' V r(r)M^f(r,r')^(r'), 
where iff = zf — zf s and M is the matrix 

8 2 S S 2 S 

8z?(t)8z'/(t') 8z9(t)Sz^(t') 
8 2 S S 2 S 

Sz*"(r)8zf(T') 8zr<*(T)5zP(T>) 



(3.14) 



(3.15) 



with elements taken at the stationary fields. 

Inserting (|3.14| in (|3.12p we obtain, after integrating 
over the z-fields, 



[Z n ] av = e ~ nSclt [ a s ' m2 • A l' <! s r: ' ( T ~ T '^ Xs< - T ~ T,S) } 
where 

S cS = S [af , m 2 , ^ , qf , Xs ) - - In Det M, 
We conclude, because of (13.3)) . that 



(3.16) 



(3.17) 



F=±S cB [af,m 2 ,Al,qf(r-r'), Xs (r~r')] . (3.18) 

In the previous expressions, the subscript s means that 
the fields are taken at their stationary values: A s (m 2 = 
2i\ s is the spin gap), A s>fl = 0, Xs(t-t') and qf p {r-r'). 

The staggered magnetization af is given in terms of 
the CP 1 fields as 

1 n 1 n 



a=l 



a=l 



From the average free-energy one can derive the phase 
diagram of the system, which is shown in Fig.l 0, [2j. 
This presents a critical line separating the SG and PM 
phases, which starts, at T = 0, in the quantum critical 
point 



A 

Po= 2^ 



1 

1 + - 

7 



1 + 1 ln(l + 7 ) 



(3.20) 
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where 



7 = 37r( Ajj =— — (3 ' 21) 



and A = 1/a is the high- momentum cutoff. For p > po, 
there is an ordered (AF) Neel phase at T = 0. 

The Edwards-Anderson (EA) order parameter, which 
is used for detecting the occurrence of a SG phase 0,13, 
is given by 



(3.22) 



where 



?—X2<l af3 ("n = 0). 
n(n — 1) 



(u> n are the Matsubara frequencies). 




0.96 0.98 



Ps / PO (0) 



FIG. 1: Phase diagram for 7 = 10 2 ,A = 10 3 . The critical 
curve valid near the QCP po (solid curve). po(0) = A/2tt is 
the QCP of the pure AF system. The value ascribed to A is 
a realistic one in K (A — * 4p-A; v s : spin- wave velocity). 



IV. QUANTUM SKYRMION CORRELATION 
FUNCTIONS 



A. The General Method 

Given a theory containing an abelian gauge field in 
two dimensions, the topological excitations are magnetic 
vortices [TH, Skyrmions, which are topological ex- 

citations of the NLSM, accordingly, appear as magnetic 
(in the ^4 M -field) vortices in the CP 1 formulation. The 
correlation functions of the corresponding vortex quan- 
tum creation operator p(x, r) are obtained by treating 
this operator as a disorder variable, dual to the order 
parameter of the system, as we saw in section II. Then, 
a method of quantization has been developed [HI, Il4j |. 
where all correlation functions of the topological excita- 
tion creation operator, can be obtained by modifying the 
integrand of the partition function by adding to the cor- 
responding field intensity tensor F^, an external partic- 
ular field configuration B tiv {z]x^y). We implement this 
method below for the present case. 

We will take the CP 1 field z, as the "order" field. From 
(|3.19|) . it is clear that whenever (zf) ^ 0, the staggered 
magnetization will not vanish as well. We can then in- 
troduce the dual algebra relating zf with the topological 
excitation creation operator, in the form of (|2.2[) . Using 
the fact that the symmetry group of the model is U(l), 
we have [l3[ 

/•*(*, *)*f(y> *) = CX P {* arg(y - x)} zf{y, t)p(x, t). 

(4-1) 

The vortex creation operator p satisfying (|4.ip is given 

by 



fi(x, t) = exp I i2n f d 2 r arg(r - x) ^ [z?*ir? 

{ J a=l 



«]M)j>, 
(4.2) 



where irf is the momentum canonically conjugate to 
zf. Whenever the A^-field kinematics is described by 
a Maxwell term, this can also be written as 



A detailed study of the phase diagram of the sys- 
tem, explicitly showing the J and AJ dependence can 
be found in 

We would like to emphasize that in obtaining the above 
phase diagram, both in the CP 1 and NLSM versions, 
quantum fluctuations have been included in the deriva- 
tion of the average free-energy, hence it goes beyond the 
mean-field approximation. The fact that the ordered AF 
phase is washed out from any T ^ is a clear evidence 
for this. 

Our extended NLSM depends on the coupling J pre- 
cisely in the same way as the usual NLSM. Thus, the 
analysis of the validity of the approximations made, as a 
function of the value of this coupling follows closely the 
corresponding analysis well-known in the simple NLSM 



At(x,t) = exp<^ z2tt / d 2 r arg(r - x)djF Ui (r, t) 



by just considering the field equation. Then, using the 
Cauchy-Riemann equation for arg(r — x) and the analyt- 
ical properties of this function, we obtain, equivalcntly 



/i(x, t) = exp < i2n I 



d^F^{r,t) 



(4.3) 



where F^ v is the field intensity tensor corresponding to 
A, 

This last form of the operator is the most useful for 
obtaining the correlation functions. It is not difficult to 
infer from (|4.3|) that, for a Maxwell- type action, for which 
S g is quadratic in F^ u , the /i-field two-point correlation 
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function will be given by [131 ] 

(fj,(x,T)(J(y,T')) = J X>yl M exp \~S g F M „ + B^ v {z] x, y) 
where 



(4.9) 



(4.4) 



B fiV (z;x,y) = 2n [ V ^ ' ^6^*^ - 0- ( 4 - 5 ) 

J a:— (x,r ) 

This turns out to be a general expression for the topolog- 
ical excitation correlators, namely, it holds for any action 
S g [Fftv] whatsoever, irrespective of its form (13, 

Now, an important adjustment must be made, in or- 
der to adapt this result to the present disordered sys- 
tem. In order to recover the physical thermodynamics, 
we must take the limit where the number of replicas van- 
ishes, n — ► 0, hence we must be careful when defining the 
physical correlator in this limit. From (|4.2p . we see that 
we actually have 



(4.6) 



because fields belonging to different replicas commute. 
It is natural, therefore, to define the physical correlation 
function as the geometrical average among the n replicas, 
before taking the limit n — > 0, namely 



(W^phys = lim 



n^O 



limW)* (4.7) 

n— >0 



.Q = l 

where the correlator on the r.h.s. is given by (|4.4p . 

B. The Effective Gauge Field Theory 

We must now determine the form of the A^-field ac- 
tion S g [F^ v ] in l|4.4|) . in order to calculate the topological 
quantum correlation function. For this purpose, we start 
from the CP 1 description of the system, [2| expand the 
action in the Zi's and fields up to the second order 
around the stationary points (zf = A^ = 0) and per- 
form the quadratic Z{ integrals. We must stress that, 
expanding around zf = 0, we are only considering the 
case when (a) = and therefore this analysis does not 
apply to the ordered AF phase. This has been already 
considered elsewhere fla. 



VzVz*VA^e 



-S[zf ,zf * ,A M ,m 2 ,Xs(r— r').?"' 3 (t-t')] 

(4.8) 

We took the fields A, x an d Q af3 at their stationary point 
leaving the integrals in z\ and A^. Then, integrating on 
the z^-fields and expanding the action in the A^-field up 
to quadratic fluctuations around the stationary configu- 
ration A^ — 0, we get 



where, 



S g [A»] = - / dr \d\A 



S 2 S cS [A, 



bAJA v 



A v . 



Au.=0 



(4.10) 

In Appendix A, we show that the action for the A^- 
field can be written as 



-- J dr' I dV^(r,r)n(r,r';r,r')F^(r',r') 



where k = fj, a = q (5^7) and 
n(fc, u n ) 



(4.11) 



T 



dx 



L \k\ 2 +io 2 ]Jo \k\ 2 x(l — x) + xlu 2 + m 2 

(4-12) 

Notice that the a-term in the effective action (|4.11l) is 
different from zero only in the SG phase, where q~o =/= 0. 



C. The Skyrmion Correlation Functions 

We may write (|4.4p as 

(v(x,T)^(y,T')) = [ VA 







j VA„ exp < 









-A^[-nn\A, 



(4.13) 



where 



Id = k + all 



(4.14) 



Introducing a gauge fixing term, we can integrate 
(|4.13p in A^, obtaining 

(fi(x , t)hH$, r'))ph ys = exp S^tt 2 jf d£, x jf d Vp e x ^ a e^ a 



x d^Ftf- m 6 - m) - \b^ [u]^ } . ( 4. i r, i 



where x = (x, r),y = (y, r'), 



F(x;r)=Tj2 



;F(k,uj n y K - x e 



ik-x —iuj n T 



(4.16) 
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and 



F(k,w n ) 



n 

^0 



aT 



dx 



(4.17) 



[\k\ 2 +ujI} 2 Jo \k\ 2 x(l -x) +xul + m 2 
Using the identity 

in (I4.15| we can see that the first term cancels the last 
one in (|4.13p (and in (|4.15[l ) and the second one gives 



B. The Paramagnetic Phase 

In the PM phase, we have both a = and q Q = and, 
consequently, only the k term of the function F in (|4. 17|) 
contributes to the skyrmion correlation function (|4.18D . 
We calculated this term, for equal times, in Appendix B. 
From (IB~6l) we have, at T = 0, 



(Mfl(^)/4(£ T ))ph^ 



PM T^O 



exp 



\x - y\ 



(5.4) 



For an arbitrary temperature, T, we obtain the large- 
distance behavior (see (|B.7|) ) 



(Ufa r)^(y, rO)p hy8 = exp {4n 2 [F(x - y; r - r') - F(e; 0)]} , 

(4.18) 

where e is a short-distance cutoff. Introducing the renor- 
malized skyrmion field operator 

[j, R (x, t) = exp {27r 2 F(e; 0)} pfi, r), 

we get the renormalized and finite correlation function 

(/i fi (f,r)/x|j(y,r'))piiy S = exp {Att 2 F(x - y; r - r')} , 

(4.19) 

In Appendix B, we calculate the inverse Fourier trans- 
forms of the k and a terms of the function F. These 
allow us to determine the large distance behavior of the 
quantum skyrmion correlation functions (|4.18|) . 



LARGE DISTANCE BEHAVIOR OF 
SKYRMION CORRELATORS 

A. Neel Phase 



\x—y\ — >oc 



exp 



\x - y\ 



ttT\x - 



y\ 



V2 



cothlv^TrTl. 



y. 



From the above expression, we may infer that 

tt 2 kT' 



(Mi 



exp 



2V2 



?0. 



(5.5) 



(5.6) 



This nonzero result is the one to be expected in the PM 
phase where the topological excitations should not be 
genuine excitations, due to the absence of spontaneous 
symmetry breaking. The above result just confirms this 
fact, by stating that the topological excitation quantum 
state \hr) pm is not orthogonal to the ground state, being 
therefore trivial. 



The skyrmion two-point correlation function has been 
evaluated in the ordered AF phase, which occurs on the 
line T — 0; p s > po in [lfj. It presents the following 
large distance behavior 



(p R (x, T)n R (y,r))^ ys 



\x—y\ — >oo 



exp {— 2tt(j 2 \x — y\} 



(5.1) 

where a is the staggered magnetization satisfying a 2 = 
g [p s — Po] 111 ■ The equation above implies 



{m) AF = 0, 

meaning that the quantum skyrmion states \/ir) af 
orthogonal to the ground state and are, consequently, 
nontrivial. The exponential decay of the skyrmion cor- 
relator, conversely, implies the corresponding quantum 
excitations (skyrmion) have a gap E g = 2ira 2 . Notice 
that this gap, which may be written as 



(5.2) 



are 



E 9 = T [P* - Po] 



(5.3) 



vanishes as we approach the quantum phase transition to 
the SG phase at the quantum critical point (T — 0, p s = 
Po)- 



C. The Spin-Glass Phase 

In the SG phase, we have a = but now qo ^ 0. Then, 
the a term of the function F in (|4.17[) will contribute to 
the skyrmion correlation function (|4.18[) .Wc have calcu- 
lated this term, for large distances and equal times, in 
Appendix B. This has a logarithmic behavior and there- 
fore dominates the large distance behavior of the function 
F. Indeed, according to (|B.9|) . we have 



F(x-y;0) 



\x— y \— *oo 



<loD 



24Tr 2 p s m 2 



1 - 



3T 2 
Ira 2 



where C is a constant. 

Inserting this result in (|4.19|) , we get 

/ /-» \ t /-* w ST 1 $\ — >oc 1 

{pR{x,T)p R {y,T)} phys > 



x-y\ 



hxC\x-y\, 
(5.7) 



(5.8) 



where 



qoD 
6p s m 2 



1 - 



3T 2 



2m 2 



(5.9) 
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In realistic systems, we always have T 2 <C m 2 (typically 
T ~ lOi-T and m ~ 100-ftT, therefore f is always positive 
in the SG phase. As a consequence of this 

<M^)/4(y,T)>^ y8 (5.10) 

and therefore 

(fi R ) SG = 0. (5.11) 

This result shows that the quantum skyrmion states are 
orthogonal to the ground state in the SG phase, being 
therefore non-trivial quantum states. The power-law be- 
havior of their two-point correlator, however implies that 
they have a zero excitation gap. This is a quite inter- 
esting result. It reveals the existence of a new dual- 
ity relation between the AF and SG phases. The usual 
order-disorder duality occurs between the AF and PM 
phases, namely the order and disorder parameters, (a) 
and (/it) are respectively nonzero in each of these two 
phases, while the other vanishes. 

The duality between the AF and SG phases, for both 
of which qEA 7^ 0, however, is of a different nature. The 
spin excitations are gapless whereas the topological ones 
are gapped in the AF phase and conversely the spin exci- 
tations are gapped whereas the topological ones are gap- 
less in the SG phase. Our finding of gapless topological 
excitations in the SG phase is a key step in establishing 
this new duality relation. 

D. The Spin-Glass Phase and the BKT Mechanism 

There is also an important fact related to the exis- 
tence of gapless topological excitations and their asso- 
ciate power-law correlators in the SG phase. When us- 
ing the CP 1 language the skyrmions become vortices. 
The power-law behavior of their correlation functions is 
a clear indication that we have a Berezinskii-Kosterlitz- 
Thouless (BKT) two-dimensional system of vortices [§[ . 
Indeed, quantum vortices do have a large-distance power- 
law decay in the low-temperature phase, which exists be- 
low the critical point in a BKT system [9( • This explains 
the existence of a SG phase at a finite temperature, as 
observed in 0, 0| : it is a BKT phase supporting gapless 
vortices. Conversely, the presence of gapless spin wave 
excitations in the AF phase, washes this phase away at 
any finite temperature through the well-known Mermin- 
Wagner mechanism [6]. This is the explanation for the 
asymmetry found between the SG and AF phases of this 
system, the former persisting at a finite T, whereas the 
other only remains at T = 0. 

The operation of the BKT mechanism in this sys- 
tem only occurs because it is possible to use a CP 1 de- 
scription, which presents gapless vortices in a SG phase. 
Thereby one can understand why we can have a stable 
SG phase in a quantum SO (3) disordered 2D Heisenberg 
system, but not in the corresponding Ising system 
The latter does not allow a CP 1 formulation with the 



corresponding vortices and therefore cannot display the 
BKT mechanism. 



E. A New Characterization of the Spin-Glass 
Phase 

The SG phase is the realization of one of the possible 
phases allowed by the dual algebra of spin and topological 
excitation operators, namely the one where both (a) = 
and (fx) = 0. This criterion - both order and disorder pa- 
rameters vanishing - can be used as a more symmetrical 
alternative for the characterization of a SG phase then 
the usual one where (a) = and qEA 0. We can state, 
equivalently that a SG phase is one where the quantum 
topological excitations are nontrivial gapless states of the 
Hilbert space. 

One can speculate whether this is a general property 
or a peculiarity of the present model. From the point 
of view of energetics, it is clear that the creation of a 
skyrmion defect out of an ordered ground state, such as 
the one we have in an AF phase, costs a finite amount of 
energy since a number of spins must be flipped, in order 
to create the quantum defect state. Conversely, creating 
such a defect on a disordered ground state, as the one we 
have in a paramagnetic phase, clearly does not change the 
state of the system and therefore costs no energy. The 
skyrmion operator actually does not not create a truly 
new state. We can consider that the topological defects 
are condensed in the ground state, thereby producing the 
disordered PM state. Skyrmions here are not genuine 
excitations. 

In a spin glass phase the ground state is also a disor- 
dered state. Therefore, when creating a skyrmion defect 
on such a state, we may conclude that, as in the PM 
phase, there will be no cost in energy, because of the 
disordered character of the ground state. However, dif- 
ferently from the PM phase, the SG ground state is a 
"frozen" disordered state. Consequently the new state 
generated by the inclusion of the skyrmion will be non- 
trivially different from the ground state, despite being 
also disordered and having zero energy cost for its cre- 
ation. The skyrmion state, hence, must be orthogonal to 
the ground state, implying the correlation function must 
vanish at large distances. This way the large distance 
behavior of the skyrmion correlation function detects the 
frozen nature of the ground state. The result is the oc- 
currence of zero energy nontrivial topological states or- 
thogonal to the ground state. The fact that they bear 
a nonzero topological charge clearly distinguishes them 
from the ground state. 

These arguments may be generalized for Ising systems, 
for instance, by replacing skyrmions by Bloch walls and 
can be applied whenever topological defects may be in- 
troduced. It seems to allow a broad characterization of 
a SG phase, for a vast class of systems, as one for which 
both the order and disorder parameters vanish, namely 
(<t) = 0, (fi) = 0, where a is the spin and fi is the topo- 
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logical excitation operators. 
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VI. CONCLUSION 

The disordered magnetic system considered in this 
work presents all the possible phases allowed by the du- 
ality relation, which exists between the staggered spin 
operator and the creation operator of quantum topologi- 
cal excitations, namely quantum skyrmions. There is an 
ordered antiferromagnetic phase with (a) =/= 0, (/i) = 0, 
a paramagnetic phase with (a) = 0, (//) ^ and a spin- 
glass phase, with (a) = 0, (fx) = 0. 

The PM and AF phases are dual to each other, in 
the sense that the order and disorder parameters, respec- 
tively (<j) and (fi) show a complementary behavior, being 
zero or not, respectively, in each of the two phases. An 
interesting duality relation, however, also exists between 
the AF and SG phases, concerning the gap of topological 
and spin excitations. In the AF phase, we have gapless 
spin excitations, namely m = 0, whereas the quantum 
topological excitations have a finite gap proportional to 
a 2 and given by (|5.3p . The SG phase, conversely, presents 



spin excitations with a non- vanishing gap, namely m ^ 
whereas the topological excitations are gapless, ac- 



cording to (|5.8p . The AF and SG phases, therefore are 
dual with respect to the gap of the spin and topologi- 
cal excitations. The paramagnetic phase presents spin 
excitations with a gap m =^ 0, but there are no gen- 
uine topological excitations in the Hilbert space, since 
the creation operator of topological excitations acting on 
the ground-state produces basically the same state. 

The fact that in the SG phase the skyrmion state is 
orthogonal to the vacuum in spite of having zero energy 
is an indication that the ground state is frozen albeit 
disordered. We, therefore arrive at an alternative char- 
acterization of a spin glass state, as one in which both the 
order and disorder parameters vanish and the quantum 
topological excitations in a magnetic systems are gapless 
but nontrivial. 

Our results indicate that the stability of the SG phase, 
where both (er) and (fx) vanish, at a finite temperature, 
is a consequence of the BKT mechanism. This does not 
apply to Ising systems, where a vortex picture of topo- 
logical excitations does not exist. 

The average free-energy of a spin-glass is expected to 
have a large amount of local minima. The problem of de- 
termining the absolute minimum is a difficult one. One 
can find results in the literature, indicating that renor- 
malization group flows may drive a replica symmetric 
solution towards a broken replica symmetry one. One 
possibility is that this would happen in a high coupling 
(J) limit. This would bring the system to a new local 
minimum where the replica symmetry would be broken. 
This is a very interesting subject for future investigation. 
It is, however beyond the scope of the present work. 
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APPENDIX A: EFFECTIVE ACTION FOR A L , 



Let us derive here the expression for the action of the 
A^-field. We will consider only the case a = 0, namely, 
the PM and SG phases. From (|3.17|) it follows, in this 
case that 



IP" = 



S 2 S eS [A fJ 



5AnSA v 



n 8 A M <5 A l 



[lnDetM[A 



wu u =o ■ 



(Al) 



By expanding in A^ and taking the derivatives be- 
fore making the limit A^ — > 0, we find in frequency- 
momentum space k^ = (w n , k) , 

n^(w n ,fc) = n{k 2 5^ - k"k v ) + Aq T{uo n ,k)8^ , (A.2) 
where 

d 2 q 1 



T(w n ,k) = T 



T 



( 27r ) 2 [\q\ 2 + m 2 ][\k - q\ 2 + u)\ + m? 



dx 



4ttJo \k\ 2 x(l - x) + xu)% + 



(A.3) 



m 



Multiplying and dividing the last term in (|A.2|) by k 2 = 
(\k\ 2 +uj 2 ) and adding a pure gauge term, the result (|4.11|) 
follows, with 



APPENDIX B: FINITE TEMPERATURE 
INVERSE FOURIER TRANSFORMS 



1. The Inverse Transform of 



We have T 1 [ pjgr ] g iven by 



d 2 k e i(k-£-iu n r) 

(2tt) 2 k 2 +iu 2 



(B.l) 



The k- integral may be easily done [16(, yielding 

T 



oo 



+2 ^ K 



(B.2) 



Integrating on the angular k- variable, we get 



where Kq is a modified Bessel function. For r = we can 
perform the sum [l6l |. The logarithmic term is canceled 
and we get, up to a constant (which will not contribute to 
the skyrmion correlation function, as we can infer from 

(SH) 



0) 



1 



47r|a; — y\ 



{1 



-1 d 



x J (fc|f|)e-^" r 



n(|fc|,o 



(B.7) 



n=l 



1 



1 



T 2 + u& 



(B.3) 



We immediately see that, for zero temperature we get 
t^o 1 



0) 



4ir\x — y\ 



(B.4) 



For finite temperatures, we obtain from (|B.3|) . at large 
distances 



Fi{x; 0) 

ttT\x — y 
+ 



\x-y\-nx> 



1 



1 



An\x - y\ { 2 
coth[V27rT|f - 



(B.5) 



We are actually interested in the large distance behavior 
of Ti(x; 0). In this regime, only the k — » will contribute 
to the fc-integral. Since n(|fc|,w„) is always regular for 
k — > 0, we can replace it by n(0,o->„). Then the o; n sum 
may be performed, for r = 0, yielding to leading order in 



^ 2 (s;0) 



1 



127rm 2 



dfcfc 



fc 2 + e 2 ' 
(B-8) 

The last integral gives [l6[ ifo(e|^|) — — mC|x|, hence 



3T 2 " 






lim / 


2m 2 





\x\ — VOO 



1 



127rm 2 



3T 2 
2^ 



\nC\x\ (B.9) 



2. The Inverse Transform of ^^~Tp 



Now we have JF" 



given by 



d 2 fc n(|fc|,^„) i( k.2-iu nT ) 



J 7 2(x;t)=T) i — — -— —c 

^ I (2tt) 2 [k 2 + w 2 ] 2 



(B.6) 



This immediately leads to (|5.7[) and to the power law 
behavior of the skyrmion correlation function in the SG 
phase, (|5.8p and (|5.9p . The constant C will not con- 
tribute to the skyrmion correlation function, as we can 
infer from (I4.18P and therefore is not important in this 
framework. 
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